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We study in this paper the parallel spin current in an antiferromagnetic semi- 
conductor thin film where we take into account the interaction between itinerant 
spins and lattice spins. The spin model is an anisotropic Heisenberg model. We use 
here the Boltzmann's equation with numerical data on cluster distribution obtai- 
ned by Monte Carlo simulations and cluster-construction algorithms. We study the 
cases of degenerate and non-degenerate semiconductors. The spin resistivity in both 
cases is shown to depend on the temperature with a broad maximum at the tran- 
sition temperature of the lattice spin system. The shape of the maximum depends 
on the spin anisotropy and on the magnetic field. It shows however no sharp peak 
in contrast to ferromagnetic materials. Our method is applied to MnTe. Comparison 
to experimental data is given. 

PACS numbers: 75.76.+j ; 05.60. Cd 

I. INTRODUCTION 

The behavior of the spin resistivity p as a function of temperature (T) has been shown 
and theoretically explained by many authors during the last 50 years. Among the ingredients 
which govern the properties of p, we can mention the scattering of the itinerant spins by the 
lattice magnons suggested by Kasuyan^, the diffusion due to impurities^, and the spin-spin 
correlation." First-principles analysis of spin-disorder resistivity of Fe and Ni has been also 
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recently performed .- 

Experiments have been performed on many magnetic materials ranging from metals to 
semiconductors. These results show that the behavior of the spin resistivity depends on 
the material : some of them show a large peak of p at the magnetic transition tempera- 
ture Tcr others show only a change of slope of p giving rise to a peak of the differen- 
tial resistivity dp/dT.—^ Very recent experiments such as those performed on ferromagne- 
tic SrRuOs thin films^^, Ru-doped induced ferromagnetic Lao.4Cao.6MnO^, antiferroma- 
gnetic e-(Mni_j.Fea;)3,25Ge^, semiconducting Pro.rCao.sMnOs thin filmsi^,superconducting 
BaFe2As2 single crystals^^, Lai_i.Sr2.Mn03^^ and Mni^xCTxTe^^ compounds show different 
forms of anomaly of the magnetic resistivity at the magnetic phase transition temperature. 

The magnetic resistivity due to the scattering of itinerant spins by localized lattice spins 
is proportional to the spin-spin correlation as proposed long-time ago by De Gennes and 
Friedel^, Fisher and Langer-, and recently by Kataoka^. They have shown that changing 
the range of spin-spin correlation changes the shape of p. In a recent work, Zarand et al.- 
have showed that in magnetic diluted semiconductors the shape of the resistivity versus T 
depends on the interaction between the itinerant spins and localized magnetic impurities 
which is characterized by a Anderson localization length (. Expressing physical quantities 
in terms of ( around impurities, they calculated p and showed that its peak height depends 
indeed on this localization length. 

In our previous work—"— we have studied the spin current in ferromagnetic thin films. 
The behavior of the spin resistivity as a function of T has been shown and explained as an 
effect of magnetic domains formed in the proximity of the phase transition point. This new 
concept has an advantage over the use of the spin-spin correlation since the distribution of 
clusters is more easily calculated using Monte Carlo simulations. Although the formation of 
spin clusters and their sizes are a consequence of spin-spin correlation, the direct access in 
numerical calculations to the structure of clusters allows us to study complicated systems 
such as thin films, systems with impurities, systems with high degree of instability etc. On 
the other hand, the correlation functions are very difficult to calculate. Moreover, as will 
be shown in this paper, the correlation function cannot be used to explain the behavior 
of the spin resistivity in antiferromagnets where very few theoretical investigations have 
been carried out. One of these is the work by Suezaki and Mori^'' which simply predicted 
that the behavior of the spin resistivity in antiferromagnets is that in ferromagnets if the 
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correlation is short-ranged. It means that correlation should be limited to " selected nearest- 
neighbors" . Such an explanation is obviously not satisfactory in particular when the sign of 
the correlation function between antiparallel spin pairs are taken into account. In a work 
with a model suitable for magnetic semiconductors, Haas has shown that the resistivity p 
in antiferromagnets is quite different from that of ferromagnets.— In particular, he found 
that while ferromagnets show a peak of p at the magnetic transition of the lattice spins, 
antiferromagnets do not have such a peak. We will demonstrate that all these effects can be 
interpreted in terms of clusters used in our model. 

In this paper, we introduce a simple model which takes into account the interaction 
between itinerant spins and localized lattice spins. This is similar to the s — d model^i. The 
lattice spins interact with each other via antiferromagnetic interactions. The model will be 
studied here by a combination of Monte Carlo simulation and the Boltzmann's equation. 
As will be discussed below, such a model corresponds to antiferromagnetic semiconductors 
such as MnTe. An application is made for this compound in the present work. 

The paper is organized as follows. In section [TTl we show and discuss our general model 
and its application to the antiferromagnetic case using the Boltzmann's equation formula- 
ted in terms of clusters. We also describe here our Monte Carlo simulations to obtain the 
distributions of sizes and number of clusters as functions of T which will be used to solve 
the Boltzmann's equation. Results on the effects of Ising-like anisotropy and magnetic field 
as well as an application to the case of MnTe is shown in section III. Concluding remarks 
are given in section IV. 

II. THEORY 

Let us recall briefly principal theoretical models for magnetic resistivity p. In magnetic 
systems, de Gennes and Friedel^ have suggested that the magnetic resistivity is proportional 
to the spin-spin correlation. As a consequence, in ferromagnetically ordered systems, p shows 
a divergence at the transition temperature Tc-, similar to the susceptibility. However, in 
order to explain the finite cusp of p experimentally observed in some experiments. Fisher 
and Langer^ suggested to take into account only short-range correlations in the de Gennes- 
Friedel's theory. Kataoka^ has followed the same line in proposing a model where he included, 
in addition to a parameter describing the correlation range, some other parameters describing 
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effects of the magnetic instability, the density of itinerant spins and the apphed magnetic 
field. 

For antiferromagnetic systems, Suezaki and Mori-9 proposed a model to explain the ano- 
malous behavior of the resistivity around the Neel temperature. They used the Kubo's 
formula for a.n s — d Hamiltonian with some approximations to connect the resistivity to 
the correlation function. However, it is not so easy to resolve the problem. Therefore, the 
form of the correlation function was just given in the molecular field approximation. They 
argued that just below the Neel temperature T/v a long-range correlation appears giving 
rise to an additional magnetic potential which causes a gap. This gap affects the electron 
density which alters the spin resistivity but does not in their approximation interfere in 
the scattering mechanism. They concluded that, under some considerations, the resistivity 
should have a peak close to the Neel point. This behavior is observed in Cr, a — Mn and 
some rare earth metals. Note however that in the approximations used by Haaa^i, there is 
no peak predicted. So the question of the existence of a peak in antiferromagnets remains 
open. 

Following Haas, we use for semiconductors the following interaction 

V = Y,J{r-Rn)s-Sn (1) 
n 

where J{f— Rn) is the exchange interaction between an itinerant spin s at r and the lattice 
spin S„ at the lattice site i?„. In practice, the sum on lattice spins S„ should be limited at 
some cut-off distance as will be discussed later. Haas supposed that V is weak enough to 
be considered as a perturbation to the lattice Hamiltonian given by Eq. ( TT5|) below. This is 
what we also suppose in the present paper. He applied his model to ferromagnetic doped 
CdCr2Se4^"— and antiferromagnetic semiconductors MnTe. Note however that the model by 
Haas as well as other existing models cannot treat the case where itinerant spins, due to the 
interaction between themselves, induce itinerant magnetic ordering such as in (Ga,Mn)As 
shown by Matsukura et al.'' Note also that both the up-spin and down-spin currents are 
present in the theory but the authors considered only the effect of the up-spin current since 
the interaction "itinerant spin" -"lattice spin" is ferromagnetic so that the down-spin current 
is very small. This theory was built in the framework of the relaxation-time approximation 
of the Boltzmann's equation under an electric field. As De Gennes and Friedel, Haas used 
here the spin-spin correlation to describe the scattering of itinerant spins by the disorder 
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of the lattice spins. As a result, the model of Haas shows a peak in the ferromagnetic case 
but no peak in the antiferromagnetic semiconductors. Experimentally, the absence of a peak 
has been observed in antiferromagnetic LaFeAsO by McGuire et al.— and in CeRhlns by 
Christianson et al.— 



A. Boltzmann's equation 

In the case of Ising spins in a ferromagnet that we studied before^, we have made a 
theory based on the cluster structure of the lattice spins. The cluster distribution was incor- 
porated in the Boltzmann's equation. The number of clusters 1] and their sizes have been 
numerically determined using the Hoshen-Kopelmann's algorithm (section III Bp . We work 
in diffusive regime with approximation of parabolic band and in an s — d model. We consi- 
der in this paper that in our range of temperature the Hall resistivity is constant (constant 
density). To work with the Born approximation we consider a weak potential of interac- 
tion between clusters of spin and conduction electrons. We suppose that the life's time of 
clusters is larger than the relaxation time. As in our previous paper— we use in this paper 
the expression of relaxation time obtained from the Boltzmann's equation in the following 
manner. We first write the Boltzmann's equation for /, the distribution function of itinerant 
electrons, in a uniform electric field E 

(^)(f^) ^ (|).». (2) 

m oe at 

where p is the equilibrium Fermi-Dirac function, k the wave vector, e and m the elec- 
tronic charge and mass, e the electron energy. We next use the following relaxation-time 
approximation 

{-^)coll = ■~(:^)' fk = fk - (3) 

where is the relaxation time. Supposing elastic collisions, i. e. k = k', and using the 
detailed balance we have 



dt (27r) 



j[wk^M'-fl)W. (4) 
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where Q is the system volume, Wk'^k the transition probabihty between k and k'. We find 
with Eq. 03]) and Eq. (jl]) the following well-known expression 

1 Q f 

= / K/,fc(l-cos^)] 

xsm9k''^dk'd9d(f), (5) 

where 9 and are the angles formed by k' with k, i. e. spherical coordinates with z axis 
parallel to k. 

We use now in Eq. (|5]) the "Fermi golden rule" for Uk,k' and we obtain 



1 n 



[[cok,k'il-cos{9))]sm{9)k'^dk'd9d(P (6a) 

Tk l^^j J 

O0k,k' = ^1^1 < k'\J{r)\k > \'6{k' - k) (6b) 

where J(r) is the exchange integral between an itinerant spin and a lattice spin which is 
given in the scattering potential, Eq. ([1]). One has 

J{r)^J{\f^-R^\) (7) 

Note that for simplicity we have supposed here that the interaction potential J(r) depends 
only on the relative distance r' = \r — Rn\, not on the direction of r — i?„. We suppose in 
the following a potential which exponentially decays with distance 

J(r) = He-^/« (8) 

where Vq expresses the magnitude of the interaction and ^ the averaged cluster size. After 
some algebra, we arrive at the following relaxation time 

1_ = 32V>vr _ 1 _ {2^kf)^ 

Tk, {2khf ^ i + mff [i + {2ikf)Y ^ ^ 

where kf is the Fermi wave vector. As noted by Haas^^, the mobility is inversely proportional 
to the susceptibility x- So, in examining our expression and in using the following expression 
X = Yl^^viOr^ where ?7(^) is the number of clusters of size ^, one sees that the first term 
of the relaxation time is proportional to the susceptibility. The other two terms are the 
corrections. 
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The mobility in the x direction is defined by 

We resolve the mobility fi^ explicitly in the following two cases 
- Degenerate semiconductors 



k 

^tWe)n - 2.(ff^i^{^ml±^^r (lib) 

where D — ^-r . We arrive at the following mobihty 



mDkf 4^2 



^-^g/^=zfl^[ ' '.7.'"' ]' (12b) 

The resistivity is then 



,4y>W 4^ (13a) 



We can check that the right-hand side has the dimension of a resistivity : 



[C?[s] 

mm]. 

Non-degenerate semiconductors 
One has in this case = exp{—/3ek) 



j:f'-M^-^f'r'^'V^/2 (14a) 

k 

V k'(df2/de)n - 2n(—?/' - (—V/'\l + ^ ^^""^^ + fl4b) 

2_^i, [aj,/ae)n-Z7:[^^) 2^(4^)2/3^/12^ h'l3 ^ m/3' ^ ^ 

ne'n' ,2m/3,i/2r 2 x 16m£2 6(8m^2)2 

m2i:)(4^2j20j:^ ^2 ffjj fj4p2 

p = - (14d) 

(T 

Where D = ^1%^ 
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Note that the formulation of our theory in terms of cluster number 7] and cluster size ^ is 
numerically very convenient. These quantities are easily calculated by Monte Carlo simu- 
lation for the Ising model. The method can be generalized to the case of Heisenberg spins 
where the calculation is more complicated as seen below. In section IIII Al we will examine 
values of parameter Vq where the Born's approximation is valid. 

B. Algorithm of Hoshen-Kopelmann and Wolff's procedure 

We use the Heisenberg spin model with an Ising-like anisotropy for an antiferromagnetic 
film of body-centered cubic (BCC) lattice of A^^ x Ny x Nz cells where there are two atoms 
per cell. The film has two symmetrical (001) surfaces, i.e. surfaces perpendicular to the z 
direction. We use the periodic boundary conditions in the xy plane and the mirror reflections 
in the z direction. The lattice Hamiltonian is written as follows 

'H = jJ2^.-S,+Aj2stS] (15) 

{id) {id) 

where Sj is the Heisenberg spin at the site i, is performed over all nearest-neighbor 

(NN) spin pairs. We assume here that all interactions including those at the two surfaces 
are identical for simplicity : J is positive (antiferromagnetic), and A an Ising-like anisotropy 
which is a positive constant. When A is zero, one has the isotropic Heisenberg model and 
when A — )■ oo, one has the Ising model. The classical Heisenberg spin model is continuous, 
so it allows the domain walls to be less abrupt and therefore softens the behavior of the 
magnetic resistance. Note that for clarity of illustration, in this section (H.B) we suppose 
only NN interaction J. In the application to MnTe shown in section III.C, the exchange 
integral is distance-dependent and we shall take into account up to the third NN interaction. 

Hereafter, the temperature is expressed in unit of J/ks, ks being the Boltzmann's 
constant. A is given in unit of J. The resistivity p is shown in atomic units. 

For the whole paper, we use = Ny = 20, and Nz = 8. The finite-size effect as well 
as surface effects are out of the scope of the present paper. Using the Hamiltonian (fT5|) . we 
equilibrate the lattice at a temperature T by the standard Monte Carlo simulation. In order 
to analyze the spin resistivity, we should know the energy landscape seen by an itinerant 
spin. The energy map of an itinerant electron in the lattice is obtained as follows : at each 
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position its energy is calculated using Eq. ([S]) within a cutoff at a distance Di = 2 in unit of 
the lattice constant a. The energy value is coded by a color as shown in Fig. [T]for the case 
A = 0.01. As seen, at very low T {T = 0.01) the energy map is periodic just as the lattice, i. 
e. no disorder. At T = 1, well below the Neel temperature Tn ~ 2.3, we observe an energy 
map which indicates the existence of many large defect clusters of high energy in the lattice. 
For T ^ T/v the lattice is completely disordered. The same is true for T = 2.5 above T/v. 

We shall now calculate the number of clusters and their sizes as a function of T in order 
to analyze the temperature-dependent behavior of the spin current. 




constant a and A = 0.01, for T = 0.01, T = 1.0, T = 2.0 and T = 2.5 (from left to right, top to 
bottom, respectively). The values of energy corresponding to different colors are given on the right. 



The scattering by clusters in the Ising case in our previous model^^ is now replaced in 
the Heisenberg spin model studied here, by a scattering due to large domain walls. Counting 
the number of clusters in the Heisenberg case requires some particular attention as seen in 
the following : 

- we equilibrate the system at T 

- we generate first bonds according to the algorithm by Wolff it consists in replacing 
the two spins where the link is verified the Wolff's probability, by their larger value 
(FigH 

- we next discretize S^, the z component of each spin, into values between —1 and 1 with 
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subtract by the lattice at T=0 
only Szfcontribution 
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Figure 2: The successive steps in the appHcation of the algorithm by Wolff to the case of Heisenberg 
spin. See text for explanation. 

a step 0.1 

- only then we can use the algorithm of Hoshen-Kopelmann to form a cluster with the 
neighboring spins of the same S^- This is how our clusters in the Heisenberg case are 
obtained. 

Note that we can define a cluster distribution by each value of Sz- We can therefore 
distinguish the amplitude of scattering : as seen below scattering is stronger for cluster 
with larger Sz- We have used the above procedure to count the number of clusters in our 
simulation of an antiferromagnetic thin film. We show in Fig. [3] the number of cluster r) 
versus T for several values of Sz- 

We have in addition determined the average size of these clusters as a function of T. The 
results are shown in Fig. |H One observes that the size and the number of clusters of any 
value of Sz change the behavior showing a maximum at the transition temperature. 

The resistivity, as mentioned above, depends indeed on the amplitude of Sz as seen in 
the expression 



m 1 



m 
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300 




T 

Figure 3: Number of clusters versus temperature for anisotropy A = 0.01 (upper) and A = 1 
(lower). The values of Sz are 1, 0.8 and 0.6 denoted by circles, squares and triangles, respectively. 
Lines are guides to the eye. 

III. RESULTS 

A. Effect of Ising-like Anisotropy 

At this stage, it is worth to return to examine some fundamental effects of Vq and A. It 
is necessary to know acceptable values of Vq imposed by the Born's approximation. To do 
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Figure 4: Average size of clusters versus temperature for anisotropy A = 0.01 (upper) and A = 1 
(lower). The values of 5^ are 1, 0.8 and 0.6 denoted by circles, squares and triangles, respectively. 
Lines are guides to the eye. 



this we must calculate the resistivity with the second order Bom's approximation. 



An 



F{e, 



^[ f t/Ve-^^-V(r) - — /" d^re-'"^-^ I d'r'e 
^ Air r 



J(r')] 



we find, with D 
1 



K =\k-k'\ = k[2{l - cose)]^/2 and J(r) = Vqc 
ri327iQm 



DV'k 



1 + 



+ 



Vie 



Tk ° ^[l + {2ikYY 3[1 + (2eA:)2]2^ [l + (2eA:)2]2^^ 12(2A;2) 



(17a) 
(17b) 



(18) 
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The first term is due to the first order of Born's approximation and the second and third 
terms to corrections from the second order. We plot p{Born2) / p{Bornl) versus T in Fig. [S] 
for different values of Vq, p{Bornl) and p{Born2) being respectively the resistivities calcu- 
lated at the first and second order. We note that the larger this ratio is, the more important 
the corrections due to the second-order become. From Fig. [5l several remarks are in order : 

- The first order of Born's approximation is valid for small values of Vq as seen in the case 
Vq = 0.01 corresponding to a few meV. In this case the resistivity does not depend on 
T. This is understandable because with such a weak coupling to the lattice, itinerant 
spins do not feel the effect of the lattice spin disordering. 

- In the case of strong Vq such as Vq = 0.05, the second-order approximation should be 
used. Interesting enough, the resistivity is strongly affected by T with a peak corres- 
ponding to the phase transition temperature of the lattice. 



Figure 5: Ratio Residue= p{B or n2) / p{B or nl) versus T for Vo=0-05 (squares, upper curve) and 
0.01 (circles, lower curve). See text for comments. 

We examine now the effect A. Figure |6] shows the variation of the sublattice magneti- 
zation and of Tjv with anisotropy A. We have obtained respectively for A = 0.01, A = 1, 
A = 1.5 and pure Ising case the following critical temperatures T/v — 2.3, 4.6, 5.6 and 6.0. 
Note that the pure Ising case has been simulated with the pure Ising Hamiltonian, not 
with Eq. f|T5|) (we cannot use A = oo). We can easily understand that not only the spin 
resistivity will follow this variation of T/v but also the change of A will fundamentally alter 
the resistivity behavior as will be seen below. 
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Figure 6: Sublattice magnetization versus T for several values of anisotropy A. From left to right 
A = 0.01, ^ = 1, yl = 1.5 and pure Ising spin. 

The results shown in Fig. [7] indicate clearly the appearance of a peak at the transition 
which diminishes with increasing anisotropy. If we look at Fig. |4] which shows the average 
size of clusters as a function of T, we observe that the size of clusters of large Sz diminishes 
with increasing A. 

We show in Fig. |8]the pure Heisenberg and Ising models. For the pure Ising model, there 
is just a shoulder around T/v with a different behavior in the paramagnetic phase : increase or 
decrease with increasing T for degenerate or non degenerate cases. It is worth to mention that 
MC simulations for the pure Ising model on the simple cubic and BCC antiferromagnets 
where interactions between itinerant spins are taken into account in addition to Eq. ([1]), 
show no peak at all^*^. These results are in agreement with the tendency observed here for 
increasing A. 

B. Effect of Magnetic Field 

We apply now a magnetic field perpendicularly to the electric field. To see the effect of 
the magnetic field it suffices to replace the distribution function by 



From this, we obtain the following equations for the contributions of up and down spins 




(19) 



mc 




(20) 
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Figure 7: Spin resistivity versus T for several anisotropy values A in antiferromagnetic BCC 
system : A = 0.01 (circles), 1 (squares), 1.5 (triangles). Upper (lower) curves : degenerate (non 
degenerate) system. 

pt - 2^ (5. - 1) — —^^ — h+Aek^} 

where 5*2 is the domain-wall spin (scattering centers) and Vq is the coefficient of the exchange 
integral between an itinerant spin and a lattice spin [see Eq. dH])]. 

Figures [9] and [10] show the resistivity for several magnetic fields. We observe a split in 
the resistivity for up and down spins which is larger for stronger field. Also, we see that the 
minority spins shows a smaller resistivity due to their smaller number. The reason is similar 
to the effect of A mentioned above and can be understood by examining Fig. [TT] where we 
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Figure 8: Spin resistivity for pure Heisenberg (circles) and Ising (squares) models in antiferroam- 
gnetic BCC system. Upper (lower) curves : degenerate (non degenerate) system. 

show the evolution of the number and the average size of clusters with the temperature in a 
magnetic field. By comparing with the zero-field results shown in Figs. [3] and HI we can see 
that while the number of clusters does not change with the applied field, the size of clusters 
is significantly bigger. It is easy to understand this situation : when we apply a magnetic 
field, the spins want to align themselves to the field so the up-spin domains become larger, 
critical fluctuations are at least partially suppressed, the transition is softened. 
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Figure 9: Resistivities of up (circles) and down (squares) spins versus T for two magnetic field's 
strengths in the degenerate case. Top (bottom) : B = 0.6(1.5). 




T 



Figure 10: Resistivities of up (circles) and down (squares) spins versus T for two magnetic field's 
strengths in the non degenerate case. Top (bottom) : B = 0.6(1.5). 

C. Application to MnTe 

We have chosen a presentation of the general model which can be applied to degenerate 
and non-degenerate semiconductors and semi metals. The application to hexagonal MnTe is 
made below with the formulae of both degenerate and non-degenerate cases, for comparison. 
Hexagonal MnTe has a big gap (1.27 eV), but it is an indirect gap. So, thermal excitations 
of electrons to the conduction band may not need to cross the gap channel. This may justify 
the use of the degenerate formulae. In the degenerate case, kf depends only on the carrier 
concentration n via the known formula '■ kj = (37r^ra)^/^ . We use for MnTe n = 2x lO^^cm"^ 
mentioned below. For the non-degenerate case, kf is not necessary. Note that in the case of 
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Figure 11: Upper : Number of clusters, Lower : Average size of clusters, versus T for several values 
of Sz and for magnetic field B = 1.5. Circles : Sz = 1, squares : Sz = 0.8, triangles : 5^ = 0.6. 
Lines are guides to the eye. 

pure intrinsic semiconductors, fcj is in the gap and its position is given by the law of mass 
action using parabolic band approximation. In doped cases, band tails created by doped 
impurities can cover more or less the gap. But this system, which is disordered by doping, 
is not a purpose of our present study. 

In semiconductors valence electrons can go from the valence band to the conduction 
band more and more as the temperature increases. Therefore, the carrier concentration is a 
function of T. Our model has a number of itinerant spins which is independent of T in each 
simulation. However in each simulation, we can take another concentration (see Ref. 19) : 
the results show that the resistivity is not strongly modified, one still has the same feature. 
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except that the stronger the concentration is the smaller the peak at Tc becomes if and only 
if interaction between itinerant spins is taken into account. Therefore, we believe that generic 
effects independent of carrier concentration will remain. Of course, the correct way is to use 
a formula to generate the carrier concentration as a function of T and to make the simulation 
with the temperature-dependent concentration taking account additional scattering due to 
interaction between itinerant spins. Unfortunately, to obtain that formula we have to use 
several approximations which involve more parameters. We will try this in a future work. 

In the case of Cdi_a;Mna;Te, the question of the crystal structure, depending on the doping 
concentration x remains open. Cdi_a.Mna;Te can have one of the following structures, the 
so-called NiAs structure or the zinc-blend one, or a mixed phase.—""— 

The pure MnTe crystallizes in either the zinc-blend structure^'' or the hexagonal NiAs 
one^ (see Fig. [T^ . MnTe is a well-studied p-type semiconductor with numerous applications 
due to its high Neel temperature. We are interested here in the case of hexagonal structure. 
For this case, the Neel temperature is T/v = 310 K— . 




Figure 12: Structure of the type NiAs is shown with Mn atoms only. This is a stacked hexagonal 
lattice. Up spins are shown by black circles, down spins by white ones. Nearest-neighbor (NN) 
bond is marked by 1, next NN bond by 2, and third NN bond by 3. 



The cell parameters are a = 4.158A and c = 6.7lA and we have an indirect band gap of 



Eg = 1.27eV. 



Magnetic properties are determined mainly by an antiferromagnetic exchange integral 
between nearest-neighbors (NN) Mn along the c axis, namely Ji/ks = —21.5 ±0.3 K, and a 
ferromagnetic exchange ^2/^5 ~ 0.67 ±0.05 between in-plane (next NN) Mn. Third NN in- 
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teraction has been also measured with J^/ks — — 2.87±0.04 K. Note that the spins are lying 
in the xy planes perpendicular to the c direction with an in-plane easy-axis anisotropy^^. 
The magnetic structure is therefore composed of ferromagnetic xy hexagonal planes anti- 
ferromagnetically stacked in the c direction. The NN distance in the c direction is therefore 
c/2 ^ 3.36 shorter than the in-plane NN distance a. 





Figure 13: Number of clusters (upper) and cluster size (lower) versus T for MnTe structure 
obtained from Monte Carlo simulations for several values of 5"^ : 1 (circles), 0.8 (squares), 0.6 
(triangles). Lines are guides to the eye. 

We have calculated the cluster distribution for the hexagonal MnTe using the exchange 
integrals taken from Ref.— and the other crystal parameters taken from the literature^"—. 
The result is shown in Fig. [T31 The spin resistivity in MnTe obtained with our theoretical 
model is presented in Fig. [T3] for a density of itinerant spins corresponding to n = 2 x 10^^ 
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cm~^, together with "normahzed" experimental data. The normahzation has been made by 
noting that the experimental resistivity R in Ref. 41 is the total one with contributions 
from impurities and phonons. However, the phonon contribution is important only at high 
T, so we can neglect it for T < 310K. While for the contribution Rq from fixed impurities, 
there are reasons to consider it as temperature-independent at low T. From these rather 
rude considerations, we extract Rq from R and compare our theoretical with R — Rq. This 
is what we called " normalized resistivity" in Fig. [HJ 




600 



T(K) 

Figure 14: Normalized spin resistivity versus T in MnTe : theoretical non-degenerate case (circles), 
theoretical degenerate case (squares) and experimental results (stars) from Chandra et al4i. Expe- 
rimental data lie on the degenerate line for T > 140 K. See text for comments. 

Several remarks are in order : 

i) the peak temperature of our theoretical model is found at 310 K corresponding the the 
experimental Neel temperature although for our fit we have used only the above-mentioned 
values of exchange integrals 

ii) our result is in agreement with experimental data obtained by Chandra et al.^^^ for 
temperatures between 140 K and 280 K above which Chandra et al. did not unfortunately 
measured 

iii) at temperatures lower than 140 K, the experimental curve increases with decreasing T. 
Note that many experimental data on various materials show this 'universal' feature : we can 
mention the data by Li et al}^, Du et al.— , Zhang et al.^^, McGuire et al.^^ among others. Our 
theoretical model based on the scattering by defect clusters cannot account for this behavior 
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because there are no defects at very low T. Direct Monte Carlo simulation shows however that 
the freezing indeed occurs at low T both in ferromagnets^i^ and antiferromagnets^ giving 
rise to an increase of the spin resistivity with decreasing T. There are several explanations 
for this experimental behavior among which we can mention the fact that in semiconductors 
the carrier concentration increases as T increases, giving rise to an increase of the spin 
current, namely a decrease of the resistivity, with increasing T in the low-T region. Another 
origin of the increase of p as T — )■ is the possibility that the itinerant electrons may be 
frozen (crystallized) due to their interactions with localized spins and between themselves, 
giving rise to low mobility. On the hypothesis of frozen electrons, there is a reference on the 
charge-ordering at low T in Pro.5Cao.5Mn05r^ due to some strain interaction. A magnetic 
field can make this ordering melted giving rise to a depressed resistivity. Our present model 
does not correspond to this compound but we believe that the concept is similar. For the 
system Pro.sCao.sMnOa, which shows a commensurate charge order, the "melting" fields at 
low temperatures are high, on the order of 25 Tesla^^. 

iv) the existence of the peak at T/v = 310 K of the theoretical spin resistivity shown in 
Fig. [13] is in agreement with experimental data recently published by Li et al.^^ (see the 
inset of their Fig. 5). Unfortunately, we could not renormalize the resistivity values of Li et 
al.— to put in the same figure with our result for a quantitative comparison. Other data on 
various materials^'^'^^ also show a large peak at the magnetic transition temperature. 

To close this section, let us note that it is also possible, with some precaution, to apply our 
model on other families of antiferromagnetic semiconductors like CeRhlns and LaFeAsO. 
An example of supplementary difficulties but exciting subject encountered in the latter 
compound is that there are two transitions in a small temperature region : a magnetic 
transition at 145 K and a tetragonal-orthorhombic crystallographic phase transition at 160 
K.— 1^ An application to ferromagnetic semiconductors of the n-type CdCr2Se4^ is under 
way. 

IV. CONCLUSION 

We have shown in this paper the behavior of the magnetic resistivity p as a function 
of temperature in antiferromagnetic semiconductors. The main interaction which governs 
the resistivity behavior is the interaction between itinerant spins and the lattice spins. Our 
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analysis, based on the Boltzmann's equation which uses the temperature-dependent clus- 
ter distribution obtained by MC simulation. Our result is in agreement with the theory 
by Haas^i : we observe a broad maximum of p in the temperature region of the magnetic 
transition without a sharp peak observed in ferromagnetic materials. We have studied the 
two cases, degenerate and non-degenerate semiconductors. The non-degenerate case shows 
a maximum which is more pronounced than that of the degenerate case. We would like 
to emphasize that the shape of the maximum and its existence depend on several physical 
parameters such as interactions between different kinds of spins, the spin model, the crystal 
structure etc. In this paper we applied our theoretical model in the antiferromagnetic semi- 
conductor MnTe. We found a good agreement with experimental data near the transition 
region. We note however that our model using the cluster distribution cannot be applied at 
very low T where the spin resistivity in experiments is dominated by effects other than s — d 
scattering model of the present paper. One of these possible effects is the carrier proliferation 
with increasing temperatures in semiconductors which makes the resistivity decrease with 
increasing T experimentally observed in magnetic semiconductors at low T. 
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